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Abstract. We present follow-up spectroscopy of emission line candidates detected on near-infrared narrow band 
images in the inner Galaxy (IHomeier et al. 2003L The filters are optimized for the detection of Wolf-Rayet stars 
and other objects which exhibit emission-lines in the 2 /im region. Approximately three square degrees along the 
Galactic plane have been analyzed in seven narrow-filters (four emission-lines and three continuum). We have 
discovered 4 new Wolf-Rayet stars and present coordinates, finding charts, and K-band spectra. 
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1. Introduction 

Dust in the plane of our Galaxy severely compromises any 
attempt at complete samples of young stars, but can be 
mediated by observing at longer wavelengths. In Homeier 
et al. (2003), we discussed the difficulty in conducting 
searches for evolved massive stars in the optical passbands, 
and detailed our observing strategy for finding obscured 
emission-line stars using narrow band filters near 2 /jm. 
The eventual goal of such a project is the identification 
of this young stellar population to increase the number of 
known Wolf-Rayet (WR) stars, to improve the statistics of 
WR stars in clusters vs. isolated (likely runaways) WRs, 
and to use WR stars as a tracer of Galactic structure. In 
this paper we report on spectroscopic follow-up from our 
imaging program in the inner Milky Way. 
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2. Observations & Reduction 

For details of our reduction strategy for the survey images 
and candidate selection we refer the reader to Homeier et 
al. (2003). Our candidates are selected with a combination 
of photometric index and image subtraction. We avoided 
the region within 5' of the Galactic center due to severe 
crowding. 

In Table 1 we present central wavelengths 
and FWHMs for our chosen set of K-band filters 
i|Blum fc Damineli 19991 IHomeier et al. 2003|l . Four fil- 
ters are centered on the characteristic WR wind emission 
lines of He I 2.06 //m, C IV 2.08 /im, H I Br 7 2.166 ^m, 
and He II 2.189 /im, and the additional three continuum 
filters are at 2.03 /xm, 2.14 /xm, and 2.248 /im. Thus each 
line filter measurement has a continuum point to the red 
and blue, which accounts for spatially variable extinction. 

Spectra were taken with SOFI l|Moorwood et al. 1998| l 
on the European Southern Observatory's (ESO) 3.6m New 
Technology Telescope (NTT) at La Silla on the nights of 
2002 May 1 - 2,4 - 6 and with OSIRIS 1 on the Cerro 
Tololo Inter- American Observatory's (CTIO) Blanco 4m 
during the nights of 2000 July 12-13 and 2002 July 9-12. 
Seeing varied between 0.8" and 2" during the run at the 

1 OSIRIS (Ohio State Infrared Imager and Spectrograph) 
is a collaborative project between Ohio State University and 
CTIO. Osiris was developed through NSF grants AST 90-16112 
and AST 92-18449. OSIRIS is described in the instrument man- 
uals found on the CTIO Web site at http:/ /www. ctio.noao.edu 
See also DePoy et al. 1993. 
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Table 1. Filter Description 



Central A (/im) 


FWHM (/im) 


2.032 


0.010 


2.062 


0.010 


2.077 


0.015 


2.142 


0.020 


2.161 


0.022 


2.191 


0.013 


2.248 


0.024 



NTT, between 1 - 1.5" during the 2000 run at the Blanco 
4m, and between 0.8" and 1.5" during the 2002 run at the 
Blanco 4m. 

SOFI and OSIRIS deliver plate scales of 0.29" pixel" 1 
and 0.40" pixel -1 , respectively. With SOFI we used the 
GR Grism Red for a wavelength coverage from 1.5— 2.4/im, 
and the the //3 spectroscopic mode on OSIRIS for a wave- 
length coverage of 2.0 — 2.4/im. The spectral resolution at 
2.2 m is X/SX w 800 for SOFI and w 1200 for OSIRIS. 

We offset our targets by 5" along the slit and sub- 
tracted these paired observations from each other to re- 
move the background. Spectra were extracted with the 
APALL task in IRAF 2 . AO - 2V stars were observed at 
a range of airmass, the Br7 line was removed by lin- 
ear interpolation, and telluric correction of the science 
spectra was performed with these spectra. For the NTT 
observations, wavelength calibration was performed us- 
ing arc spectra. For the 2000 Blanco observations, wave- 
length calibration was accomplished using OH emission 
lines ( |01iva fc Origlia 1992| ) . For the 2002 Blanco observa- 
tions, rough wavelength calibration was accomplished us- 
ing telluric features in the standard star spectrum. Finally, 
the science spectra were normalized to the continuum. 

3. Results: 4 New Wolf-Rayet Stars 

We have discovered four Wolf-Rayet stars: three WC and 
one WN. Names were given to the stars following the 
convention in the Seventh Catalog of Wolf-Rayet Stars 
Ijvan der Hucht 200 l|l . In Table we list coordinates for 
the four stars and in Figures 1 1181 we present spectra and 
finding charts. 

3.1. Spectra 

The three WC stars, WRIOlp, WRIOlq, and WR102ca, 
have nearly identical 2 — 2.4/im spectra. The most signif- 
icant differences are: (1) WRIOlp lacks the C IV line at 
2.139 /im, (2) WRIOlq has a C III line at - 2.11/mi which 
is stronger than the blend of C IV lines at 2.08/im, in con- 
trast to WRIOlp and WR102ca, (3) WRIOlq has a line 
at 2.22 /im of an unknown transition, and (4) WR102ca 
appears to have more dust than WRIOlp and WRIOlq, as 

2 IRAF is distributed by the National Optical Astronomy 
Observatories 



Table 2. 



Star 


RA (J2000) 


Dec (J2000) 


WRIOlp 


17:45:42.47 


-28:52:53.3 


WRIOlq 


17:45:57.78 


-28:54:46.1 


WR102ca 


17:46:13.04 


-28:49:25.4 


WR102ka 


17:46:18.12 


-29:01:36.5 




Fig. 1. Finding chart for WRIOlp. The star H5 from 
Cotera et al. 1999 is ~ 51" to the Southwest. North is 
up and east is to the left. The image size is 5 x 5'. 
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Fig. 2. Spectrum of WRIOlp. 

its lines are weaker overall due to excess continuum emis- 
sion by dust ; see the discussion below on excess emission. 

The spectra are similar to the spectrum of star 309 
in the Quintuplet cluster, which Figer et al. (1999) clas- 
sified as earlier than WC8 due to its excess emission 
from He II 3.09 /im and lack of strong emission from He I 
2.058 /im . Based on their similar spectral morphology in 
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Fig. 3. Finding chart for WRIOlq. North is up and east 
is to the left. The image size is 5 x 5'. 
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Fig. 4. Spectrum of WRIOlq. 

the 2 — 2.4/im region, these spectra are consistent with 
the same classification, with varying amounts of dust di- 
lution. Of the published spectra of known WRs other than 
the Quintuplet WRs, the closest morphological match is 
the WC9 star WR 88 l|Eenens et al. 1991|) . We therefore 
classify these stars as WC8-9. 

The 2 — 2.4/im spectrum of WR102ka shows strong 
lines of He I 2.06/im and Bi'7, weak lines at 2.103/im 
CTII/NIII, 2.114/rai Hel/NIII, 2.189 /im Hell, and pos- 
sibly an unidentified feature at 2.139 /im. In Figure El 
we present the 2.05 — 2.25/im region of the spectrum of 
WR102ka to illustrate these weak features. 

Of note is the strength of He I 2.06/im to 
Br7 ; the HcI/Bi-7 ratio is ~ 2. A few of the 
Galactic Center He I stars have similar Hel/Br7 ratios 
H Allen et al. 19901 IKrahhe et al. 199ll IHIum et al. 19951 
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Fig. 5. Finding chart for WR102ca. The concentration of 
stars is the Quintuplet cluster, and WR102ca is a member. 
North is up and east is to the left. The image size is 5 x 5'. 
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Fig. 6. Spectrum of WR102ca. 



|Tamblyn et al. 1996| |Najarro et al. 1997| ) , but most also 
have stronger 2.11/xm features than WR102ka. The 
most similar spectral morphology is displayed by two 
LMC stars, BAT99 45 (HDE 269582, BE294) and 
BAT99 13 (Hen S9, S k -66 40) < |McGregor et al. 19881 
|Breysacher et al. 1999| ). These stars were originally clas- 
sified as Ofpe/WN9, but later re-classified as WNlOh and 
WN10 fsee ICrowther fc Smith 19971 for an explanation). 
The general picture is that stars between 40 — 60 M© 
evolve through an LBV stage to become WN8 stars, 
and in their dormant stage are identified as WN9 — 
11 stars IjCrowther fc Smith 19971 |Pasquali et al. 1999| ). 
The star BAT99 45 is at least spectrally variable. The 
Hel/Br7 strength changed by a factor of 10 between 
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Fig. 7. Finding chart for WR102ka. North is up and east 
is to the left. The image size is 5 x 5'. 
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Fig. 8. Spectrum of WR102ka. 
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Fig. 9. Spectrum of WR102ka on a different scale. 
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Fig. 10. J — H vs.i? — Ks plot with magnitudes taken 
from the 2MASS catalog; see text. The main sequence and 
giant tracks are indicated by the solid line and dotted lines 
l|Koornneef 1983JI and the reddening line is overplotted as 
the dashed line ( |Rieke fc Lebofsk~ 985). 



1987 January and 1993 September HMcGregor et al. 1988| 
IBlum et al. 1995jl . Due to the similarity of WR102ka's 
2-2.4/xm spectra with BAT99 45 and BAT99 13, we clas- 
sify it as WN10. 

3.2. JHKs Photometry 

In Table we present JHKs magnitudes from the latest 
calibration of the Two Micron All-Sky Survey (2MASS) 
catalog 3 A J— H vs. H — Ks plot is presented in FigurelTUI 
with the main sequence and giant tracks as solid and 
dotted lines IjKoornneef 1983(1 . and the reddening line for 
the bluest stars as a dashed line according to a Rieke & 
Lebofsky (1985) extinction law. Objects falling to the right 
of the reddening line have colors inconsistent with extinc- 
tion along the line of sight for the bluest stars, i.e. they 
have excess emission at Ks, usually attributable to hot 
dust. 

J, H, and Ks errors were added in quadrature to 
produce the error bars in Figure El For WRIOlq and 
WR102ca, no J errors were listed in the 2MASS catalog, 
and we have therefore not shown a J — H error. 

In FigurelTDwe plot J — H vs. H — Ks values for most 
of the known Wolf-Rayet stars in the Galaxy taken from 
the 2MASS catalog. The J, H, Ks magnitudes were re- 
trieved using the positions given in the Seventh Catalog of 
Wolf-Rayet Stars l|van der Hucht 20 01). Magnitudes are 
included if the positions in the 2MASS catalog were within 
3" of the input coordinates. If more than one object was 
found, the nearest object was taken; misidentifications 
are possible. We also plot J — H and H — K' for the 
WN and WC stars in the Quintuplet cluster, and the 
four stars with featureless if— band spectra referred to as 
"cocoon stars" , or "dusty, late-type WC stars" (DWCLs; 



http://www.ipac.caltech.edu/2mass/ 



Homeier et al.: Results from a NIR search for Galactic WRs 



5 




Fig. 11. J — H vs.H — Ks plot for a sample of Galactic WN and WC stars, the stars discovered here, and the Quintuplet 
WN, WC, and DWCL stars. The J,H,Ks magnitudes were retrieved from the 2MASS catalog using the positions 
given in the Seventh Catalog of Wolf-Rayet stars Ijvan der Hucht 2 001). Magnitudes are included if the positions in 
the 2MASS catalog were within 3" of the input coordinates. If more than one object was found, the closest was 
taken; misidentifications may be possible. J, H, and K' magnitudes for the Quintuplet stars were taken from Figer 
et al. (1999); where two measurements were made, only the first is plotted. The main sequence and giant tracks are 
indicated by the solid line and dotted lines IjKoorn neef 1983J) and the reddening line is overplotted as the dashed line 
flRieke fe Lebofsky 1985| >. 

Table 3. 2MASS Magnitudes 



Name J error H error Ks error J — H H — Ks 

WRIOlp 1(132 (U32 PXoB O059 VL20 O049 3^27 L85 - 

WRIOlq 16.26 null 13.47 0.070 11.49 0.056 2.79 1.98 

WR102ca 14.56 null 12.79 0.071 10.40 0.075 1.77 2.39 

WR102ka 12.98 0.022 10.27 0.025 8.84 0.031 2.71 1.43 



|Figer et al. 19 99). For reference, the K : Ks, and K' fil- 
ters have slightly different wavelength coverages: the Ks 
filter transmission profile starts ~ 0.04 /im bluer and ends 
~ 0.06 /im bluer than K; K' starts ~ 0.04 /im bluer and 
ends ~ 0.02 /im bluer than Ks. 

The WR stars are offset from the reddening line for 
the bluest main sequence stars towards larger H — Ks. 
To further illustrate this, we include an inset of this re- 
gion in Figure Both the WN and WC subtypes follow 



sequences parallel to the plotted reddening line due to in- 
tervening interstellar dust, with the WN stars offset by 
0.1 — 0.2 magnitudes and the WC stars by another 0.4 
magnitudes. This can be understood as the effect of the 
contributions from free- free emission in their ionized winds 
and a blackbody spectrum from hot dust around most WC 
stars (see Cohen et al. 1975 and references therein), and 
possibly line emission. WR stars have strong emission lines 
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in their NIR spectra and this may affect their broad-band 
colors. 

We have also labeled the positions of WR98a and 
WR104 in Figure im both which are associated with 'pin- 
weel' nebulae HTuthill et al. 19991 IMonnier et al. 1999|l . 
The known pinwheel nebulae arise from the interaction 
of a WC wind with that of a massive companion. The po- 
sitions of these stars along the WC sequence support the 
idea that the near- infrared colors are driven by hot dust, 
and that most WC stars in the Galaxy produce dust (al- 
though not all, as some fall in the region occupied by the 
WN stars). 

WRIOlp (WC8-9) and WR102ka (WN10) have near- 
infrared colors more typical of Galactic WN stars and 
most of the Quintuplet WR stars, as does WRIOlq (WC8- 
9), although its J — H color is uncertain due to an un- 
qualified error in its J magnitude. WR102ca (WC8-9) 
appears to have a large infrared excess and at the same 
time appears to suffer from less extinction than the rest 
of the Quintuplet WR stars. We note that it is likely that 
the J magnitude is severely in error, which could give it 
an anomalously blue J — H color and an overestimated 
K— band excess. In the 2MASS H and Ks images there 
are two stars within ~ 0.08" radius of WR102ca, while in 
the J image only one of the three stars (not WR102ca) 
is easily visible. All three stars have listings of "null" for 
their estimated J error, and thus their J magnitudes is 
likely inaccurate. 

4. Null Detections 

In this section we briefly discuss the candidates which 
turned out not to be emission-line stars upon follow-up 
spectroscopy. These stars were readily identified as giants 
by their strong CO absorption lines. Only the 2.06 and 
2.08 /im line indices were used successfully to identify WR 
stars. It is unclear if this is due to a lack of early- type WNs 
which dominate at the lines of 2.17 and 2.19 /im, a sensitiv- 
ity issue due to their typically weaker lines, or a detection 
bias. Our 2.17 and 2.19 /xm line indices had considerably 
more scatter than the 2.06 and 2.08 /im counterparts, lead- 
ing us to conclude that a detection bias is highly likely. 
We do however note that we were able to use the 2.17/im 
filter images to identify several objects with nebular Br7 
emission. 

During the first set of follow-up observations, objects 
were selected purely on the photometric indices. Due to 
the large scatter in these indices, the success rate was low 
(3 emission-line objects out of 45 candidates, including 1 
WR star). Other problems included bad pixels and resid- 
ual images from bright stars. However, this was largely 
eliminated with the image subtraction step. With the in- 
clusion of this step, we were able to identify 6 emission-line 
objects, including 3 WR stars, out of 38 candidates. We 
found that selection as a photometric candidate and a pos- 
itive identification as a candidate after image subtraction 
were separately necessary but not sufficient conditions for 
the identification of an emission- line star. 



Of the candidates which were not emission-line objects, 
most were multiple objects when observed at higher angu- 
lar resolution while positioning the slit. None of the mem- 
bers of these multiple objects were emission-line stars. All 
4 of the identified WRs were apparently isolated at higher 
angular resolution. This is a clear case for the need for ob- 
servations at higher angular resolution in the inner Galaxy 
where crowding is the dominant source of photometric er- 
ror. 

5. Discussion 

5.1. WC stars 

All three of the WC stars detected have similar late-type 
WC spectra (WC8-9) with dilution by dust. WR102ca 
is in the Quintuplet cluster, which has a total of 5 
members (including this new one) with similar spec- 
tra. It also contains 4 stars with featureless K-band 
spectra, the so-called "cocoon" stars, fFiger et al. 19 99), 
which are likely also to be late-type WCs with sufficient 
dust that the continuum emission overwhelms the lines 
IJChiar fc Tielens 2001inVWti et al. 200111. 

Dust around WC stars arises in regions of high 
density, which can be achieved in the wind-wind in- 
teraction zone with a massive companion. This has 
been observationally confirmed; there are known peri- 
odic WC dust-producers where the dust is produced 
at periastron ^Williams et al. 19871 I Veen et al. 1998bjl . 
and recently, the persistent dust-producers WR104 and 
WR98a have been shown to have pinwheel-shaped neb- 
ula HTuthill et al. 19991 IMonnier et al. 1999|) produced by 
a short-period (~ 1 yr) WR-0 binary system. It is not 
known if single WC stars can also produce dust solely 
through their clumped winds (see lVeen et al. 1998al for a 
discussion). The hint of copius dusty WCs in the inner 
Galaxy where the metallicity is likely > Z Q may be a 
clue to required conditions for dust formation around WC 
stars. 

5.2. WRs in the inner Galaxy 

In Homeier et al. (2003), an estimate for the number of 
WRs expected in the inner Galactic sample was given, 
following from a simple smooth model for the distribution 
of stars and gas in the Milky Way. They concluded that 
~ 28 WR stars should be detected within the area cov- 
ered by the imaging observations in the inner Milky Way. 
The model used for the estimate given in Homeier et al. 
(2003) was normalized to the local number density of WR 
stars assumed, and also depends on the scale length of 
the Galactic disk, and the distribution of WR absolute K 
magnitudes. 

A census of WR stars in the region covered 
by our images (but outside the Galactic center) 
shows that there are ~ 35 known WRs in the 
Arches cluster ( Nagata et al. 19951 ICotera et al. 1 996 
IBlum et al.lQ0"l"l|Figer et al. 2002D , 10 in the Quintuplet 
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cluster (Fig er et al. 1999| ), and at least 2 apparently iso- 
lated WR stars near Sgr A IjCotera et al. 1999}) . for a total 
of 47. Most of these are brighter than K= 12, which was 
the limit considered in the estimate. Already we see that 
the number of known WRs exceeds the prediction by this 
very simple model by nearly a factor of two. This suggests 
a higher star formation rate in the inner Galaxy compared 
to the local neighborhood. Near-infrared observations to- 
ward the inner Galaxy have established that this is indeed 
the case; see for example the review of the Galactic Center 
by Morris & Serabyn (1996). 

From this calculation, we might expect that no new 
WRs would be found brighter than K= 12, whereas we 
have identified 4 in this region. The addition of these 4 
stars is not numerically significant, but 3 of the 4 are no- 
table for their apparently isolated location. Either these 
stars are isolated and are thus runaways, or they are the 
most massive object of an association which has not yet 
been discovered. We find no evidence in our images for 
associations or clusters around these stars, but given our 
poor angular resolution, we cannot rule out their exis- 
tence. Also, small clusters may not live very long near the 
Galactic center. Numerical simulations predict that small 
(< 2 x 10 4 Mq) clusters will disperse sufficiently within a 
few Myr so that their surface densities will be indisinguish- 
able from that of their surroundings i|Kim et al. 2 000 
IPortegies Zwart et al. 20021 ). 

6. Conclusions 

We have discovered 4 new WR stars in the inner Galaxy: 
three WC8 — 9 stars with similar spectra and one WN10 
star. One of the WC stars is a member of the Quintuplet 
cluster, while the others are not members of any previously 
known clusters or associations, nor did we find evidence 
for any previously unknown clusters or associations. We 
conclude that narrow-band observations at K— band wave- 
lengths are a potentially efficient way to identify evolved 
massive stars in the plane of the Milky Way, and further 
higher angular resolution work is highly desirable. 
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